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We fiave studied detailed nucleosynthesis in the shoclced surface layers of an 
Oxygen-Neon-Magnesium core collapse supernova with an eye to determining 
if the conditions are suitable for r— process nucleosynthesis. We find no such 
conditions in an unmodified model, but do find overproduction of N=50 nuclei 
(previously seen in early neutron-rich neutrino winds) in amounts that, if ejected, 
would pose serious problems for galactic chemical evolution. 

Subject headings: supernovae, nucleosynthesis 



1. INTRODUCTION 



The site of the r— process has been the most enduring mystery in nucle o synthesis theory 



since the publication of the seminal papers in this field (jCameron 1 119571 : iBurbidge et al. 



19571 ). Of particular promise (in their times and for some even today) have been the many 
efforts suggesting Type II supernovae as the site with the relevant conditions arising either 



in or near the exploding core (initial 
on aspects of neutrino interactions 


y championed by 


J^H), with recent atten 


tion 


"ocused 
), or in 


Woos 


lev & Hoffman 


1992 


; Wooslev et al. 


1994 


the outer lavers (Thielemann et. al. 


1979; 


Cowan et. al. 1 


1982; 


Epstein. Colgate. & Haxton 


1988). For reviews, see Cowan. Thielemann. & Truran (1 


1991) 


Wallerstein et. al. Jl997). 
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The requisite conditions for r— process nucleosynthesis in explosive scenarios with mate- 
rial freezing out from nuclear statistical equilibrium (NSE) have been derived with the general 
understanding that particular combinations of three parameters, the entropy, electron frac- 
tion ( Yp), and expansion time scale, give rise to specific features of the solar r— abundance pat- 
tern (iQian fc Woosley Ill996l ). As the wind evolves these parameters must sweep out a range 
of conditions that produce the many features of the solar r— process abundances, particularly 
the relative heights of the 2nd & 3rd peaks (IWoosley et al. Ill994j ). Recent models of core 
collapse still fall short of the necessary conditions to explain all t he abundance features of the 



solar r— process, especially the high mass {A > 130) component (IHoffman. Woosley. &: Qian 
1997h . 



A site that has received recent interest involves explosions of 0-Ne-Mg cores in 8 



IOMq supern ovae. Successful SN explosions 



much debate (IWheeler. Cowan, fc Hi 



(iHille 



1985 



jrant. Nomoto. fc Wolff 



198 



I 



lebrandt 



N^omoto 



Baron. Cooperstein. fc Kahana 



rom such systems have been the subject of 



19981) . some efforts fou nd prompt explosions 



1984J ). others did not (IBurrows fc Lattimer 



19871 ) , with the former now being completely ruled out 



with the advent of modern treatm ents of neutrino transpo rt. Appeals to late-time neutrino 



heating have also been suggested ( Mayle fc Wilson 1988h. It shoul d be noted that nearly 



every attempt has used a common progenitor model (INomoto Ill984j ). 



Recent efforts to revive the idea of a low-F ^, low-entropy scenario for r— processing have 
inclu ded detailed nucleosynthesis calculations (IWanajo et. al. II2005I : iNing. Qian. &: Meyer 
20071 ). In the former case, an unmodified SN model (computed without neutrino transport) 
provides a very low explosion energy {Eexp ~ 0.02 B), with modestly neutron rich conditions 
{Ye,min = 0.45) and no r— process. To obtain it they artificially increased the shock heating 
term to obtain larger explosions (~ 1 B) and lower electron fractions (0.14 < Ye^mm < 0.36) 
and suggest material experiencing these conditions must be ejected to explain the principle 
r— process features. 



By contrast, iNing. Qian. &: Meyer I (120071 ) propose that the necessary conditions arise in 
the shocked C-Q layers above the Q-Ne-Mg core at the location of a very steep density gra- 
dient near the edge of the mass cut in these compact stars. This allows for rapid shock wave 
passage accompanied by high peak temperatures, giving rise to a short expansion timescale. 
Such a combination has been suggeste d as a viable r— process scenario even at modest en- 
tropy with a neutron excess near zero (j Jordan fc Meyer I |2004J ) . As with the neutrino wind 
scenario, this production would be primary. Ning et. al. appeal to this scenario to help 
bolster observational suggestions that r — process nuclei with A > 130 are preferentially as - 
sociated with a low mass SN component (llshimaru et. al. II2004J : iQian &: Wasserburg 1120071 ). 



Although the initial conditions of iNing. Qian. fc Meyer I (120071 ) were taken from a stellar 
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evolution model (iNomoto Ill984l . 119871 ). the evolution of temperature and density in the mass- 
shell trajectories was derived in a semi-analytical fashion assuming a shock speed in the regio n 



of interest approaching 10 cm s . Using shock jump relations (iMatzner &: Mckee Ill999l ) 



they derived the entropy and expansion timescale and obtained resulting nucleosynthesis that 
did exhibit characteristics reminiscent of r— process abundances. The authors conclude by 
imploring further research to test their assertions in modern finely zoned supernova explosion 
simulations of progenitors in the same mass range. 



Recently, iKitaura. Janka. fc Hillebrandt I (120061) have calculated supernovae e xplosions 



in th is mass range, considering an 8.8 progenitor star with an 0-Ne-Mg core (INomoto 



1984 ) and using a sophisticated treatment of neutrino transport. The nucleosynthesis studies 
presented here are based on an update of these (spherically symmetric) simulations for a 
revised progenitor model, in which the outer layers of the helium shell and a dilute hydrogen 
envelope were added (K. Nomoto, private communication). The new simulations also include 
an improved treatment of elect ron captures and inelastic neutrino scattering by nuclei in NSE 
(ILanganke et. al. I l2003l . 120071). The principle resu lts of these simulations will be discussed 



explosion dynamics can be found in 



in a separate publication ( Miiller Sz Janka I 2008 ). Som e relevant information about the 



Janka et. al. 



(120071 ). 



2. Nucleosynthesis in an exploding 8.8 Star 

Our results survey explosive nucleosynthesis in 32 zones extracted from an 8.8 M0 SN 
model starting at the edge of the 0-Ne-Mg core and extending though the C-0 and He 
layers. The mass-cut that defined our inner most zone had an initial (pre-coUapse) radius 
of 7.717 X 10^ cm with an enclosed mass of 1.363 Mq . The last zone studied had a radius 
at the onset of core-collapse of 1.1306 x 10^ cm with an enclosed mass of 1.376 M© . The 
amount of processed ejecta is 0.013 M0 . Exterior to this was a hydrogen envelope (70% 
H, 30% He) whose outermost zone at core-collapse was at radius 6.414 x 10^^ cm with an 
enclosed mass of 2.626 Mq . The amount to total ejecta is the difference, 1.263 Mq . 

During shock wave passage, the temperature increased dramatically and in most of 
the zones exceeded Tg = 7. As such photo-disintegration will disassemble any C or O 
into a— particles and nucleons. Each nucleosynthesis calculation started at a point in the 
expansion when the temperature had declined to Tg = 9.0 (or its maximum value if it never 
achieved this), with the starting values of density and Ye also taken from the SN model, and 
followed until the temperature declined to Tg = 0.1 (or the minimum value achieved in the 
SN model trajectory, in a few cases as high as Tg = 0.5). Initial compositions for zones 
with an initial Ye < 0.5 were assumed to be a-particles and free neutrons, otherwise free 
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neutrons and protons. Since nucleosynthesis beyond the iron group requires particle capture 
on seed nuclei that here must be built up by NSE, we define a simple expansion timescale 
(in seconds) as the time for the temperature to decline from Tg = 5.0 (or its highest value 
attained) to 1/e of this value. See Tabled] for the initial conditions. 



Our reaction networks are described in lPruet et. al. I (120051 ). Neu trino induced reactions 
for free nu cleons and a— p articles were included in all calculations (IMcLaughlin fc Fuller 
19951 . Il996l ). As reported by lNing. Qian. fc Meyer I (120071 ) . the inclusion of neutrino reactions 
was not crucial to the final results of what is here essentially explosive nucleosynthesis. 



2.1. Results 

Figure [1] shows the mass weighted production factors normalized to the total amount of 
eject a given by 



In this equation, the sum is over all 32 trajectories, Xj{i) is the ma ss fraction of nu clide 
i in the jth trajectory, j is the mass fraction of nuclide i in the sun (ILodders II2003I ). Mj 



is the mass of the jth trajectory, and M^^ = 1.263 Mq is the total mass ejected in the SN 
explosion whose energy was ~ 0.1 B (~ 0.2 B in a corresponding 2D simulation). 

We see no evidence of an r— process for the starting entropy and Yg conditions given 
in the model of Miiller & Janka (2007) principally because the expansion timescales of 
the trajectories containing the most mass are long (~ 0.1 s) and t he entropies are small 



(^ 10 — 20, see Table [I]). By contrast the same quantities derived by lNing. Qian. fc Meyer 



(120071 ) were entropy = 139 fcf,/nucleon, and an expansion ti mescale to e— fold frona Tg = 5.0 



of 0.013 sec. The explosion energy agrees with model Q3 of IWanajo et. al. I (120051 ) . but the 



Ye range is more neutron rich in the zones that dominated the nucleosynthesis. 

We do note production of species above the iron group, but they terminate for nuclei 
in the N=50 closed shell. Attractive combinations of expansion timescale and entropy (for 
the given Y^) do appear to be achieved for the trajectories farthest out, but these exhibit a 
rapidly diminishing drop in the peak shock temperature and density (the last two trajectories 
are in the extended hydrogen envelope, hence they have Ye > 0.5, they never reach NSE). In 
every calculation the free nucleons were all consumed in the build up of intermediate mass 
nuclei during the alpha- rich freeze out. See Table [T] for the final a— particle mass fraction in 
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each zone and the largest mass weighted production factor. 



Indeed, the nucleosynthesis here is dominated by trajectories 7 to 12 which all exhibit 
fairly long expansion timescales, entropies near 18, and initial values of < 0.48. Between 
them they constitute 40% of the total ejecta studied. The production of the dominant N=50 
isotopes (^°Zr in particular) is most pronounced in zones 9 and 10. All of the co-produced 
species (those made within a factor of ~ 4 of the largest production factor) are made as 
themselves with the dominant reaction flows occurring in the valley of stability. 

This result is reminiscent of nucleosynthesis seen at early times in the neutrino-driven 



wind of more traditional SNII (IWitti. Janka. fc Takahashi Ill993l : iHoward et. al. Ill993l : IWoosley et al. 
19941 ). where an overproduction of N=50 isotopes presented problems for a model that at late 
times produced a high entropy solution favorable to the r— process. The problem here may 
not appear so striking, in that ^°Zr is also made with ^^Se, but both are made at a very high 
production factor. Both ™Ge and ^^Zr anchor two quasi-equilibrium clusters near the iron 
group and N=50, their abundances determine the production of '''^Se and ^^Mo, respectively. 
The other two light p— nuclei in between, ^^Kr and ^^Sr, are not members of these quasi- 



equil ibrium clusters and owe their production to nuclear flows between them (IHoffman et al. 
19961 ). The dominant species (by mass fraction) for trajectory seven {Ye = 0.47, which made 
the most ^^Mo) are ^^Zn and ^°'^^Ni (all at 20 % by mass fraction), and ^°Zr, ^^Zn, and 



70 



Ge (all ~ 5%). As here, iHoffman et al. I (119961 ) also saw limited production of the heavier 
p— nuclei ^^Mo, ^^Ru and ^®Ru, which are made in equilibrium with ^°Zr but always at a 
level smaller than that for ^^Mo. The N=50 nucleosynthesis can be eliminated in preference 
to light p— nucl e i if Yp is constrained to be in a narrow range (~ 0.485). In the SN model of 
Woosley et al. I (119941 ). the N=50 overproduction problem persisted if the amount of material 
that experienced Y^ < 0.47 was greater than 1O~^M0 . Here it is even worse (see below). 

Interestingly, the recently developed i/p— process does co-produce the full range of light 
p— nuclei from '''^Se to ^°^Pd, especially ^^Mo, ^^Ru and ^^Ru, in quantities that could for 
the first time expla i n their solar abundan ces in early proton-rich neutrino wind models 



JPruet et. al. I bood : IProhhch et."an bood ). but they fail to co-produce ^^Mo. Details of 



the nuclear physics uncertainties affecting the re action flows that determine the s olar ratio 
for ^^Mo and ^^Mo have recently been explored (IFisker. Hoffman, fc Pruet 1120071 ). Within 
current uncertainties in the proton-separation energies of ^^~^^Rh, they suggest this ratio 
could be achieved along with co-production of all the light p— nuclei between Sr-Pd in the 
z/p-process if ^p(^^Rh) = 1.64 ±0.1 MeV. 
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2.2. Discussion 

If ejected, the values of the largest normalized production factors in these zones (~ 
1000) would cause serious problems for Galactic Chemical evolution. Considering the largest 
production factor (^°Zr), the enrichment by a single event relative to solar would be 



Zr 



log 



Yzr/M„ 



(2) 



where Yzr 



- single V 

(Zr) 

Y^- Xj{Zr) X Mj is the mass yield of ^^Zi (1.75 x 1O~^M0 for the sum of zones 



7-12 where ^°Zr is produced with a mass fraction grater than 10 Mmix is th e mass of ISM 



into which the ejecta of this single event mixes (~ 300 0M^. iThornton et. al. I ( 19981 )). and 
Xq is the mass fraction of ^°Zr in the sun (1.53 x 10~^, iLodders I (120031 )). This is equivalent 
to the definition of the production factor (Eq. [1]) with M*^j replaced by Mmix- For ^°Zr, the 
enrichment is 3.8, i.e. nearly four times solar, implying a very rare event. 

If Mmix is typical of the mass that can be contaminated by such a SN and assuming 
one occurred at most once at some average location over the age of the Galaxy then the 
frequency fsN for it to occur is given from 



fsNMmixtCal 



expected 



M, 



(3) 



gas 



where Mgas ~ 10^°Mq is the mass of the gas in the Galaxy, and tcai ~ 10^° yr its age, giving 
a frequency of once every three thousand years. This is in disagreement with recent best 
estimates that suggest ex plosions of 0-Ne-Mg cores comprise 4% of all core-collapse events 
(jPoelarends et. al. 1120071 ) . which translates (assuming a present day upper hmit of two SNII 
per century) to 0.0008 yr~^. If one appeals to their firm upper limit of 20% (based on the 
many uncertainties in modeling the progenitors), then this rate could rise to 0.004 yr~^. 

Over the history of the Galaxy, this enrichment is 



Zr 



log 



all 



YzrfsNtcal 



M, 



/XeiZr) 



gas 



(4) 



which gives 10-50 times the solar value. We note that both Mgas and fsN inay have been 
different in the past, but the ratio of the two which is the rate of SN per unit mass of gas 
may vary much less as the SN progenitor formation rate is proportional to the gas available. 
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Conclusions 



We have studied in detail the nucleosynthesis in the shocked outer layers of an 0-Ne-Mg 
Type II SN to test th e asser tion that this might be a viable site for the r— process. While 
Ning. Qian. &: Mever I (120071 ) derived shock conditions from the same progenitor model we 
use ( iNomoto I (Il984j ) plus the recent supplementation with a H-envelope) in a semi-analytic 
manner with an assum ed shock velocity), our conditions arise from a detailed simulation 
(jMiiller fc Janka 1 120081) . We find cond itions for nucleosynthesis that are very different from 



those of iNing. Qian. fc Meyer I (120071 ). especially in regards to peak temperatures and ex- 
pansion timescales. Consequently we cannot support their assertion that this is a potential 
site for the r— process as we see neither the requisite conditions nor products. Our results 
are however very rem iniscent of previous calculations studying the nucleosynthesis in the 
neutrino driven wind (IWoosley et al. I Il994j : iHoffman et al. I Il996l ) . 



The combination of low entropies (~ 20) and electron fractions {Yf. < 0.47) in 5.5 x 
10^'^ Mq of ejecta from our model suggest that this one event would produce nearly four 
times the solar abundance of ^*^Zr. This is unlikely, as observations of Zr are well established 
in metal poor stars where enrichment to this level has not been seen. If this represents a 
very rare event, it leaves little room for the production of Zr from other sources (such as 
the neutrino winds of ordinary SNII), and still less for the s— process. But the initial mass 
function suggests there are many stars with masses bet ween 8-10 M^;^, and recen t studies 
indicate that 4% of these should end their lives as SNII (jPoelarends et. al. 1120071 ). We are 
therefore forced to consider that this is not a typical event (or even a very rare one) and 
that some aspect of the progenitor model, the explosion model, or the nuclear physics used 
in the determination of the nucleosynthesis is grossly in error. 

On the later point, since the reaction flows producing the abundant species move along 
the valley of stability, we feel confident that the ingredients that went into the calculation 



of the nuclear reaction cr oss sections and the particle sep aration energies (so important in 
QSE) are on firm ground ( Rauscher fc Thielemann 1120001 ). There is no complication due to 



neutrino nucleosynthesis as the radii of the zones considered, lOO's of km, are too large for 
a substantial neutrino fiuence. 

On the stellar modeling side our explosion model is base d on a very modern treatment 
(iKitaura. Janka. &: Hillebrandt I l2006l : iMiiller fc Janka 1 120081 ) , with issues of fallback being 
negligible. This leaves the pre-SN model, which was originally calculated as a He-core 
with a H-envelope subsequently added to it. There are numerous difficulties in accurately 
calculating the pre-SN structure, includ ing thermal pulses of th e unstable He shell source, 
mass loss, and dredge- up to name a few (jPoelarends et. al. 1120071 ). We view this as the most 
pertinent area where improvements should be made in understanding this important class of 
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stars and consider such improvements necessary to address the assertion that they can serve 
as crucibles for half the species above iron. 

We appreciate many valuable conversations with Yong-Zhong Qian and S. E. Woosley. 
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Fig. 1. — Mass weighted production factors characterizing the integrated nucleosynthesis 
in all 32 zones studied. Nuclides of a given element are connected by solid lines, a diamond 
surrounding the data point indicates production chiefly by a radioactive progenitor. The 
horizontal lines represent a band of co-production for nuclei made within a factor of two and 
four of ^°Zr. 



Table 1. Outflow Characteristics 



(s) Mass'' X(q)'= '^Zp^^^ 



1 


1.08 


9.2 


30 


.53 


.046 


5.00 


.74 


45 Sc 


2.95 


2 


1.09 


9.5 


29 


.52 


.051 


5.00 


.72 


45 Sc 


1.98 


3 


1.23 


9.0 


27 


.51 


.058 


10.00 


.67 




3.44 


4 


1.32 


9.1 


25 


.50 


.064 


10.00 


.59 


62 Ni 


2.04 


5 


1.47 


9.1 


24 


.50 


.073 


10.00 


.55 


62 Ni 


7.05 


6 


2.14 


9.0 


21 


.48 


.092 


5.00 


.38 


64 Zn 


45.53 
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2.27 


9.1 


20 


.47 


.100 


5.00 


.28 


■^4Se 


1001.58 
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2.40 


9.0 


19 


.47 


.107 
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.17 


■^4Sc 
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9 


2.46 
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18 


.46 
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10 
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16 
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«Se 


629.45 


12 


2.79 


9.0 


14 


.47 
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.092 


10.00 


.23 
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4.68 
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9.0 


10 


.50 
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6.00 


.24 


60Ni 


.82 


17 


2.32 


9.0 


10 


.50 


.038 


1.00 


.31 


63Cu 


.29 


18 


2.37 


9.1 


13 


.50 


.032 


1.00 


.43 


63Cu 


.43 


19 


2.66 


8.0 


14 


.50 


.025 


1.00 


.49 


63 Cu 


.49 


20 


2.90 


7.1 


15 


.50 


.017 


.10 


.58 


63 Cu 


.06 


21 


2.98 


7.0 


16 


.50 


.015 


.10 


.60 


63 Cu 


.06 


22 


3.00 


6.6 


16 


.50 


.014 


.05 


.62 


63Cu 


.03 


23 


3.07 


6.3 


17 


.50 


.014 


.05 


.64 


63Cu 


.03 


24 


3.09 


5.8 


20 


.50 


.013 


.05 


.40 


63Cu 


.01 


25 


3.32 


4.9 


42 


.50 


.013 


.05 


.00 


B2Cr 


.00 


26 


3.74 


4.3 


53 


.50 


.013 


.05 


.00 


40Ca 


.01 


27 


4.77 


3.2 


88 


.50 


.013 


.01 


.02 


«4Sr 


.05 


28 


5.24 


2.8 


108 


.50 


.013 


.01 


.00 


41 K 


.00 


29 


5.99 


2.5 


132 


.50 


.013 


.01 


.80 


39 K 


.00 


30 


6.50 


2.3 


146 


.50 


.012 


.01 


1.00 




.00 


31 


7.21 


2.0 


185 


.84 


.010 


.01 


.33 


51 V 


.00 


32 


9.07 


1.4 


381 


.85 


.006 


.01 


.30 


2 H 


.00 



''Initial conditions at radius ry (in 10^ cm) when T9 ~ 9.0 (or its peaJj value). 

''Mass of the zone in units of 10^4j\^g. The mass interior to traj. 1 was 1.363 Mq . 

•^Final a— particle mass fraction, the nucleus with the largest mass weighted production 
fa<;tor, and its value. 



